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ABSTRACT 


Clearcutting of lodgepole pine stimulates production of 
understory vegetation that may provide a grazing resource 
for livestock and big game for an estimated 20 years or more. 
Peak production of 800 to 1,000 lbs. /acre occurred about 11 
years after clearcutting. Prediction equations were developed 
to permit gross estimates of potential production and its 
change with time. Although palatability of the indigenous 
understory vegetation is low, the large acreages of lodgepole 
pine that are harvested annually warrant efforts to improve 
the quantity and quality of forage by reseeding. 


Introduction 


Lodgepole pine (Pinus contorta Dougl.) occupies 5.4 million acres in Montana (Wil- 
son and Spencer 1967). This vast resource received limited commercial use until the 
early 1950's; since then, cutting for poles, pulp, and lumber has increased to an 
estimated 7,000 acres annually. ! 


Herbaceous growth develops after clearcutting and provides a potentially important 
grazing resource for approximately 20 years while the new timber stand is developing. 
Livestock and game animals now graze some clearcut areas, but the grazing potential 
appears to be considerably underutilized. This incomplete use is the result of several 
administrative factors, chief among which are a limited knowledge of the grazing re- 
source, and a tendency to ignore it because of its transitory nature. 


Efficient management of these ranges depends largely upon our ability to predict 
the expected composition, yield, and persistence of forage on clearcuts from selected 
environmental factors. Findings from a study aimed at providing the basis for predic- 
tion are presented here. 


Methods 
FIELD AND LABORATORY METHODS 


This study was limited to unthinned, lower elevation (below 7,500 feet) clearcut- 
and-burned areas in Montana, east of the Continental Divide. Twenty-five clearcut 
sites within the Lewis and Clark and the Gallatin National Forests were selected torep- 
resent a range in post-clearcut ages (2 to 17 years had elapsed since cutting) and a 
variety of aspects. 


On each of these sites, the percentage of aerial crown cover was estimated for 
each understory plant species within thirty circular plots which were 4.8 sq. ft. in 
size and distributed in a random pattern over a l-acre macroplot. All herbaceous vege- 
tation within the circular plots was clipped at ground level near the end of the grow- 
ing season, dried at 60° C. for 24 hours, and weighed. The center of each circular 
plot served also as the center of a milacre plot for data on tree regeneration. 


Physical and chemical characteristics of the A and B soil horizons were determined 
from samples taken from two soil pits on each of the 25 sites. 


We approximated site index for lodgepole pine for 19 of the 25 sampled sites from 


measures of site index, as described by Alexander (1966), on immediately adjacent, un- 
cut stands that closely matched the clearcut areas in soils and topographic features. 


1From files of USDA Forest Service, Northern Region, Missoula, Montana. 


ANALYSIS 


Successive, full-screen,* multiple regression sets enabled us to identify combina- 
tions of environmental factors that substantially affect forage production. Final 
acceptance of such combinations was conditioned by agreement of the algebraic signs of 
estimated effects with those of expected? effects, under least squares fit. In these 
analyses, average dry weight in pounds (or yield) per acre was the dependent variable, 
and 37 environmental factors constituted the complete set of independent variables 
(listed below): 

A. Time (Number of Years Since Cleareutting) 

B. Tree Factors (Regeneration): 

Density of Lodgepole Pine, determined by actual count 
Density of All Trees, determined by actual count 
Percent Canopy Cover of Lodgepole Pine 


Percent Canopy Cover of All Trees 


Stocking Rate of Lodgepole Pine, expressed as a percent of milacre plots 
in which a tree occurs 


Stocking Rate of All Tree Species 


Average Height of Lodgepole Pine Reproduction, estimated to nearest 
1/2 foot 


Average Height of All Trees 
C. Stte Index 
D. Topographic Factors: 

Percent Slope 

Elevation 

Cosine of Azimuth +1, expressed in 2 7 radians 

Cosine of Azimuth +sine azimuth +2, expressed in 2 7 radians 
E. Sotl Factors: 


Each of the following represents two variables, one for the A horizon and one 
for the B horizon: 


Field Capacity, estimated by moisture retention at 1/3-atmosphere of pres- 
sure, using standard pressure membrane apparatus and procedures outlined 
by Richards (1954). 


*The dependent variable is fitted as a function of each of all combinations of 
specified groups of independent variables via computer. R*, the Coefficient of Deter- 
mination, is printed out for each regression. 

3"Expected" from prior knowledge. 


Wilting Point, estimated by moisture retention at 15-atmospheres of pressure, 
using standard pressure plate apparatus. 


Percent Sand, determined by hydrometer method. 
Percent Silt, determined by hydrometer method. 
Percent Clay, determined by hydrometer method. 
Available Potassium, determined by flame photometer. 


Percent Organic Matter, measured by Walkely-Black wet digestion chromic acid 
method (Jackson 1958). 


Total Nitrogen, measured by the micro-Kjeldahl method, using CuSO, and Na,SO, 
digestion mixture. 


Bulk Density, measured by the paraffin coating technique. 
Soil Reaction, or pH, measured with a glass electrode in a saturated paste. 
Horizon Thickness, measured in the field. 


Water Holding Capacity, estimated by the difference between field capacity 
and wilting point, and weighted by horizon thickness, rock content, and 
percent fine material. 


Results and Discussion 


Understory yields ranged from about 100 to 1,200 lbs./acre (dry weight), and 
averaged 644 1bs./acre. 


TIME 


Understory production is expected to be low, both immediately after logging and 
after canopy closure of the reproducing lodgepole pine stand. Between these two points 
in time, production should peak. This trend was strongly evident within the 17-year 
age span of the plots, and was approximated by the sigmoid function, e-kK, where 
e = 2.718 and k = 0.01667 (age of cut minus age at peak production)*. So the time 
effect was retained as a component of all models examined. 


Total understory production reached its peak about 11 years after clearcutting. 
The several vegetal classes, however, performed differently with time (fig. 1). Sedge 
was the only class that showed no discernible pattern; its mean is shown in figure l. 
The regression for grasses was not significant (P> .05), but the data suggest slight 
peaking at about 13 years. The regressions for forbs and shrubs were significant 
(P< .05). Forb production peaked at 10 years. Shrub production increased over the 
entire range of time, but the increase was slower after 11 years than it was pre- 
viously, and there is a slight but unconfirmed hint of a flattening of production at 
about 15 or 16 years after clearcutting. These production patterns were determined on 
Sites supporting unthinned lodgepole pine reproduction. Observations in other clearcuts 
where reproduction had been thinned suggest that thinning may delay the peaking of 
understory production a few years. If this is true, a grazing potential of perhaps 20 
years, or more, may be expected on clearcuts where tree regeneration is thinned. 


900 


TOTAL—> 


600 


300 


PRODUCTION (Lbs./acre) 


YEARS SINCE CLEARCUTTING 


Figure 1.--Trends in understory production on lodgepole pine clearcuts as 
affected by time. 


The concurrent increase in shrub production and decrease in herbaceous plants is 
attributed to grouse whortleberry (Vaccintwn scoparium) which constituted about one-half 
of the shrub cover on these sites. The roots of grouse whortleberry lie mostly within 
2 or 3 inches of the soil surface. This soil zone is freer of roots of the other vege- 
tation than are the deeper zones; also, it is often the sole recipient of moisture from 
light summer rains. The competitive advantage thus accorded grouse whortleberry enables 
it to prosper at the expense of the other understory plants. 


Except for scattered individuals of other shrub species, grouse whortleberry 
is often the sole shrub in older stands of lodgepole, where it often dominates the 
understory. 


TREE FACTORS 


Tree regeneration variables are highly unstable quantities for a few years after 
clearcutting; thus, they are of limited usefulness for predicting yields of understory 
vegetation. They are of interest, however, because they do affect these yields. 


The general trend of production with time conforms to expectation. As tree cano- 
pies thicken and root systems expand, understory vegetation succumbs to the increasing 
competition for light and moisture, as evidenced by the usual paucity of understory 
plants in overstocked stands of lodgepole pine that are at the sapling-or-older stage 
of development. 
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In this study, lodgepole regeneration reached the early sapling stage on only 
three clearcuts, the youngest of which was 14 years old. Herbage yields, however, 
started to decline on 10- and 1ll-year-old clearcuts where tree heights averaged only 
2 feet. Apparently, then, lodgepole pine adversely affects herbaceous plants while still 
in the seedling stage. Moreover, we believe competition for light was negligible on 
the sampled sites because the maximum tree canopy cover was only 37 percent and averaged 
a mere 7 percent. After discounting the time effect, tree canopy and understory yield 
were not significantly correlated. 


Very likely the roots of seedling trees had developed sufficiently by 10 years of 
age to gain a competitive advantage over herbaceous plants for available moisture. 
Horton (1958) reported that open-grown lodgepole exhibited about a 1:1 ratio of maximum 
lateral root length to stem height for the first 20 years. If the 1:1 ratio applies 
to trees in this study, in which stem heights averaged 2 feet, roots should have radi- 
ated as far as 2 feet from the stem. Because lodgepole pine density averaged 5,800 stems 
per acre, then uniformly distributed trees would be spaced about 33 inches apart. Thus, 
the roots of a given tree would overlap those of its immediate neighbors by about 1 
foot. Apparently, then, we could apply to our study sites Horton's statement that "in 
a fully-stocked condition, all of the available rooting space is presumably utilized so 
that a complex interweaving network of roots occurs."' Also, we could assume that this 
network of tree roots would compete effectively with understory plants for soil moisture. 


Of course, trees were not uniformly distributed on the clearcut sites, but the 
added effect of stocking rates on accountable variation in forage yields was suffi- 
ciently high to suggest considerable overlapping of root systems: 


R2 
yield = f (time) 0.38 
yield = f (time + % stocking of lodgepole) ~47 
yield = f (time + % stocking of all trees) .55 


Further additions of tree density or of tree canopy cover, or both, failed to raise 
R2 above 0.56. 


SITE INDEX 


Ideally, a prediction model should be based on variables for which measurements 
are easily obtained. Accordingly, we examined the suitability of site index for pre- 
dicting understory yield. For this prediction, we used data from the 19 clearcuts for 
which site index could be approximated from the immediately adjacent uncut stands. 


After fitting the time effect, site index accounted for 1 percent of the variation 
in yield. As an additive or interacting term with additional variables in the model, 
the contribution of site index to accountable variation was still too small to warrant 
its retention in prediction models. Apparently this measure of site potential for tree 
growth was not closely related to the production of associated understory vegetation. 


TOPOGRAPHIC FACTORS 


As sole additive components to the "yield = f (time)" equation, slope accounted 
for an additional 20 percent of the variation, elevation for 11 percent, and aspect for 
liess than“) “percent. 


This response to slope is interesting. One usually finds increasing slope steep- 
ness associated with shallower and rockier soils and greater runoff, and therefore with 
decreasing herbage production; but in this study, contrary to expectations, understory 
production responded positively to slope. In addition, the data revealed a strong 
coincidence of the higher slope values with the more moist sites whose expected influence 
on yield is positive. This confounding of slope with other variables cast doubt on the 
validity of using slope for predictive purposes in other areas; thus, it was eliminated 
from further consideration. 


The remaining topographic variables, elevation and aspect, were added to the "yield 
= f (time)" equation and this resulted in R2 values as follows: 


yield = f (time) 0.38 
yield = f (time + elevation) .49 
yield = f (time + elevation + aspect) Sl 


When added to more complex models, elevation did not contribute sufficiently to R? 
values to warrant its retention. 


The effect of aspect was tested two ways: as cosine azimuth + 1, which places 
minimum expected yields at due south; and as cosine azimuth + sine azimuth + 2, which 
places minimum expected yields at southwest. Discounting the time effect, neither of 
these components contributed materially (less than 2 percent) to accountable variation. 


SOIL FACTORS 


The influence of moisture on understory production was emphasized by the additive 
contribution of field capacity of the A and B horizons to the "yield = f (time)" equa- 
tion, with R2 = 0.74. The further addition of sand content of the A horizon, another 
moisture-related variable, produced the equation 


G=208 
Yield = -251.4000 + 383.5379 e 60 + 3.1646FC, 


Sand,— 55> 


+ 16.8020FC, + 291.7550 e 1,200 


with R2 = 0.79. Contributions of other soil variables added to this model were 
negligible. 


Wilting point could be substituted for field capacity in this and in alternative 
models with little sacrifice in R* values. However, the addition of wilting point did 
not improve any model containing field capacity. 


Available potassium (K) in the soil solum ranked second only to the field capacity 
and ye Tene point variables in degree of correlation with understory yields, and pro- 
duced an R* of 0.73 when added to the "yield = ff (time)'' equation. The further addi- 
tions of sand content of the A and B horizons raised the R* to 0.81. However, equally 
satisfactory results were obtained with the potassium and sand values of the B horizon 
only: 


R2 
yield = f (time) 0.38 
yield = f (time + Ky) .70 
yield = f.(timer+ Ka a Sand, ) . 80 
as in the equation 
s (i=11)" 
Yaielidy=)—=228-051 + 588.614 e 60 #702 Ke 


2.5 
and), -55 


#°282.785 € 200 


Holmes and Tackle (1962) found exchangeable potassium in the B horizons to be 
noncorrelated with heights of dominant lodgepole pine trees (r = 0.055). In this study, 
potassium in the B horizon was not significantly correlated with site index for lodge- 
pole pine (r = 0.37), but potassium in the A horizon was significantly correlated 
(P< .05) with site index (r = 0.46). A modest correlation with site index plus a 
marked effect on understory yields suggests that fertilizing with potassium may be 
beneficial both to understory production and to growth rates of lodgepole pine. 


The two above equations may serve as interim predictors of understory yield. The 
response surfaces of these equations over the full ranges of time and sand content are 
shown in figures 2 and 3 for the extremes of field capacity and of potassium, respec- 
tively. In both instances, the change between surfaces is linear between the two 
extremes shown for the third independent variable. Being additive, the models may not 
fit too well near the left edge of the surfaces, since the expectation would be for 
these to come closer together near time 0. Nonetheless, they do emphasize the substan- 
tial effect of field capacity and of potassium on understory yields. 
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Figure 2.--Response surfaces at extremes of field capactty, showtng the relations 
between understory productton, time after clearcutting, and percent sand. 
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Figure 3.--Response surfaces at extremes of potassium content, showing the relations 
between understory production, time after clearcutting, and percent sand. 


FORAGE VALUE 


Major species of the understory vegetation are given in table 1. We computed 
species composition on each site from percentage aerial crown cover. Separate palata- 
bility ratings for cattle, sheep, and elk were assigned to each of the species that 
contributed 3 percent, or more, to the cover on at least one site. Ratings ranged from 
0 to 6, with a 0 rating for unpalatable plants and a 6 rating for excellent species. 
Species contributing less than 3 percent of the cover were arbitrarily rated 3, about 
midway in the scale. Sources consulted for information on palatability included Dayton 
(1960), Hayes and Garrison (1960), Hermann (1966), and USDA Forest Service (1937). 


Table 1.--Composition, aerial crown cover, and frequency of vegetatton on 
lodgepole pine clearcuts. Includes only those species contrib- 
uting 3 percent or more to the vegetal cover on at least one site 
and occurring on at least 4 of the 25 sampled sttes 


Species :~ No.of: % Frequency : % Composition :% Aerial cover 
SaceSa Avg.: Range.: Avg. : Range : Avg. :* Range 
FORBS 
Antennarta racemosa 16 22 3-53 2 Tr-6 1 fr=2 
Arnica cordifolta 23 64 13-100 7 1-24 94 Tr-9 
Aster spp. 22 Bi 3-100 6 Tr-23 2 Tr-1l 
Eptlobtum angusttfoltum 22 47 3-100 9 Tr-26 5 Tr-9 
Erigeron spp. 5 19 7-37 2 Tr-6 1 Tr-2 
Erythrontum grandiflorun 9 41 7-97 1 Tr-5 tr Tr-1 
Fragarta spp. 20 22 3-67 3 Tr-12 1 Tr-5 
Galtun boreale 7 20 3-80 1 Tr-6 Tr Tr-3 
Hteractun cynoglossotdes 24 42 3-83 Z Tr-7 Tr Tr-2 
Heuchera cylindrica 7 10 3-40 1 Tr-5 Tr Tr-1 
Linnaea borealis 5 29 3-60 5 Tr-8 2 Tr-4 
Lupinus spp. 24 47 5-97 TOY et Tr=S8 4 Tr-20 
Taraxacum officinale 20 17 3-53 1 f te ai ig Tr-2 
Thalictrum ruderale 5 16 3-43 3 cra i Tr-3 


GRASSES AND SEDGE 


Carex geyert 25 63 37-93 10 2-42 3 Tr-9 
Agrostts spp. 14 24 3-53 3 Tr-14 1 Tr-11 
Calamagrostts rubescens 14 69 7-100 15 Tr-36 5 Tr-12 
Calamagrostis spp. 7 15 3-70 5 TR-16 2 Tr-7 
Phleum pratense 4 22 3-50 2 Tr-4 Tr Tr-2 
Trisetum sptcatun 16 34 3-90 2 Tr-5 1 Tr-2 
TREES AND SHRUBS 
Aretostaphylos uva-urst 7 19 3-50 4 Tr-17 1 Tr-4 
Juntperus communis 10 10 3-23 3 Tr-10 1 Tr-4 
Rosa spp. 12 29 3-70 2 Tr-5 1 Tr-2 
Sptraea betulifolia 23 52 7-87 vy 1-17 Z Tr-7 
Vaecetntum scopartum 22 74 30-100 18 2-36 6 Tr-17 


The percent composition of each understory species was multiplied by its palatabil- 
ity rating, and the sum of these resultant products was divided by 100 to arrive at an 
average palatability rating for all understory vegetation. Based on the 0 to 6 scale, 
average palatability ratings were approximately 1.7 for cattle, 2.0 for sheep, and 2.4 
for elk--in other words, fair for cattle and sheep, and only slightly better for elk. 
Because of the widespread occurrence and appreciable amounts of lupine on clearcut 
areas, sheep grazing might be hazardous. 


Conclusions 


Clearcutting, and thinning of subsequent lodgepole pine regeneration, stimulates 
forage production that may provide a grazing resource for livestock and big game for an 
estimated 20 years or more. A harvest of 10,000 acres/year of lodgepole pine in Montana 
appears probable in the near future. Thus, the lodgepole pine type represents a graz- 
ing potential of nearly 1/4 million acres on a sustained basis. This impressive grazing 
potential warrants an upgrading of the less impressive yields (average peak production 
was about 800-1,000 lbs./acre) and palatability (fair) of the indigenous plant cover. 
For purposes of such upgrading, it is desirable to conduct seeding studies on ways to 
increase forage yields and palatability without adversely affecting timber management 
objectives. 


Where indigenous vegetation is to be grazed, gross estimates of potential produc- 
tion, and its change over a period of time, are possible with the interim prediction 
equations presented here. These estimates permit an appraisal of grazing potentials in 
advance of timber cutting and probably will allow better integration of grazing use on 
clearcuts with that on permanent range. 


If seeding of clearcuts is contemplated, the prediction equations will be help- 


ful in the judicious allocation of monies to those sites offering maximum returns on 
investment. 
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Headquarters for the Intermountain Forest and 
Range Experiment Station are in Ogden, Utah. 
Field Research Work Units are maintained in: 


Boise, Idaho 

Bozeman, Montana (in cooperation with 
Montana State University) 

Logan, Utah (in cooperation with Utah 
State University) 

Missoula, Montana (in cooperation with 
University of Montana) 

Moscow, Idaho (in cooperation with the 
University of Idaho) 

Provo, Utah(in cooperation with Brigham 
Young University) 
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